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Abstract. Recent observations of the Lyman-break galaxy (LBG) luminosity function (LF) 
from z ~ 6-10 show a steep decline in abundance with increasing redshift. However, the 
LF is a convolution of the mass function of dark matter halos (HMF)-which also declines 
sharply over this redshift range-and the galaxy-formation physics that maps halo mass to 
galaxy luminosity. We consider the strong observed evolution in the LF from z ~ 6-10 in 
this context and determine whether it can be explained solely by the behavior of the HMF. 
From z ~ 6-8, we find a residual change in the physics of galaxy formation corresponding 
to a ~ 0.5 dex increase in the average luminosity of a halo of fixed mass. On the other 
hand, our analysis of recent LF measurements at z ~ 10 shows that the paucity of detected 
galaxies is consistent with almost no change in the average luminosity at fixed halo mass from 
z ~ 8. The LF slope also constrains the variation about this mean such that the luminosity of 
galaxies hosted by halos of the same mass are all within about an order-of-magnitude of each 
other. We show that these results are well-described by a simple model of galaxy formation 
in which cold- flow accretion is balanced by star formation and momentum-driven outflows. If 
galaxy formation proceeds in halos with masses down to 10^ Af©, then such a model predicts 
that LBGs at z ~ 10 should be able to maintain an ionized intergalactic medium as long as 
the ratio of the clumping factor to the ionizing escape fraction is C//csc ^ 10. 



Keywords: galaxies: evolution - galaxies: formation - galaxies: high-redshift - cosmology: 
observations ~ cosmology: theory 



1 Introduction 



Surveys with the WFC3 camera aboard the Hubble Space Telescope have now collected large 
samples of galaxies at redshifts z > 6 using the Lyman-break technique [1-6]. Detections 
have even been made of J-band dropouts out to z ~ 10 [7, hereafter Bll] with additional 
constraints placed on the abundance at brighter magnitudes [8] . The UV luminosity functions 
(LFs) extracted from these data clearly evolve significantly from z ~ 6-10. For example, Bll 
found one candidate J-dropout above their detection limit in the two-year Hubble Ultra 
Deep Field 09 dataset, yet nine were expected assuming no evolution in the LF from z = 8. 
Fitting an artificial Schechter function to the observed J-dropout LF yields M^y = —18.3 
when fixing a normalization and slope similar to those found from z ~ 4-6 [9]. Corresponding 
values determined at z = 7 and 8 show a dimming M^y with increasing redshift. However, 
the extent to which these values elucidate the physics of galaxy formation is unclear. 

In the context of hierarchical structure formation, high-redshift Lyman-break galaxies 
(LBGs) are hosted by dark matter halos whose abundance drops rapidly from z ~ 6-10, 
mirroring the behavior of the observed LF. In theoretical models of the LF, the halo mass 
function (HMF) is convolved with a mapping from halo mass to galaxy luminosity that 
encodes the physics of galaxy formation. We seek to determine whether any residual change 
is required in this physics to account for the observed changes in the LF or if its behavior can 
be understood entirely through that of the HMF. Such a change is anticipated in a cold-flow 
model of galaxy formation [e.g., 10, 11] where halos at higher redshift have higher growth 
rates than the same mass halos at lower redshift and so may be expected to be more luminous. 
Fluctuations around the mean accretion rate also imply a range of galaxy luminosities for a 
given halo mass in contrast to many analytic models of the mass-to-light ratio [e.g., 12, 13]. 

While hydrodynamical simulations of high-redshift galaxies [e.g., 14, 15] can, in prin- 
ciple, answer these questions, they are more computational expensive and reliant on prede- 
termined physical relations than are analytic models. We use an analytic LF model that 
was initially developed by a companion paper [16, hereafter MLll] to consider the minimum 
mass of halos that can host galaxies. Here, however, we include the most recent LF data 
and improve upon our fitting method as we focus on the redshift evolution of the galaxy 
mass-to-light ratio. In §2, we briefly review this simple model and its free parameters. We 
describe the fitting procedure in §3 and consider the change in the average luminosity for 
fixed halo mass out to z ~ 10 in §4. §5 explores fluctuations around this mean luminos- 
ity, while §6 investigates the effects of observing along the light-cone given the determined 
redshift evolution of the mass-to-light ratio. We then compare our results with theoretical 
predictions of galaxy fueling and recent numerical studies in §7 and §8, respectively. Finally, 
in §9, we discuss implications for future observations. 

2 Modeling the luminosity function 

The MLll model for the high-redshift galaxy luminosity function has three components: (1) 
the Sheth-Tormen HMF [17], (2) a cutoff halo mass below which galaxies cannot form, and 
(3) a distribution function for the luminosity of a galaxy hosted in a halo of a given mass. 
Since the HMF is well-understood, we briefly summarize the second two parts of the model 
below and refer the reader to MLll for further details. 

We assume a redshift-invariant critical suppression mass Msupp is required for a halo 
to form a galaxy. This threshold is a free parameter that contains the physics of feedback 
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processes, such as supernovae and photoionization, that prevent gas from condensing into 
lower mass halos. However, since halos are built up from those of smaller mass, we expect 
some suppression of galaxy formation in halos more massive than Mgupp as well. Using a 
merger-tree algorithm based on the excursion-set formalism to quantify this effect, ML 11 
found that the "active fraction" eAFj i-e. the probability that a halo of a given mass Mh can 
host a galaxy, is redshift-independent and well-approximated by: 

logio eAF(Mh, Msupp) = -(2.231ogio[Mh/Msupp])"2-6. (2.1) 

The active fraction is distinct from the duty cycle, which is the fraction of an object's lifetime 
that it is active enough to be observable and controls the overall normalization of the LF. 
By contrast, the active fraction we define here, uniquely specified by the ratio of halo mass 
to the free parameter Msupp, is the fraction of halos with galaxies that have ever been active 
and describes the faint-end cutoff in the LF. For a fixed fraction of gas turned into stars per 
dynamical time, an instantaneous burst of star formation gives an exponentially declining 
SFR as a function of time. Therefore, if a halo is counted in the active fraction, if it has ever 
been forming stars, then it will ever after have non-zero SFR and luminosity. 

For these halos large enough to be actively forming stars, the luminosity emitted by the 
hosted galaxy is a log-normal distribution with mean Lc and standard deviation cl: 

where 



is proportional to halo mass. The shape of equation 2.2 differs from the one assumed by Ref. 
[18], who adopted a linear Gaussian profile, but is consistent with the conditional luminosity 
function approach elsewhere in the literature [e.g., 19, 20]. This distribution results from an 
excursion-set analysis of possible merger histories for halos of a given mass. For example, a 
galaxy that has had a recent merger (or, equivalently, a rapid influx of gas) will be brighter 
than one whose last major burst of star formation was further in the past. Our treatment 
obviates the need for a traditional duty cycle; if a halo is massive enough it will be forming 
stars. Galaxies that have used up much of their gas in star formation and are too faint to be 
seen at the time of observation are simply in the tail of their luminosity distribution. Instead 
of correcting our LFs for observability by eliminating faint galaxies from the calculated LF, 
as has been typical in previous studies with an occupancy or duty cycle approach [e.g., 18], 
these sources simply contribute to fainter luminosity bins beyond the detection threshold. 
This interpretation is consistent with the results of numerical simulations [15]. Thus, our 
method replaces the common degeneracy between mass-to-light ratio and duty with a reliance 
on a"L, which describes the mixing of galaxies with different halo masses among luminosity 
bins. In this picture of the luminosity distribution function, eAF is effectively the fraction 
of halos lying in a second peak to the distribution with zero luminosity. Assuming a simple 
model for the star formation in a disk with a constant fraction of gas turned into stars per 
dynamical time, a Monte Carlo analysis of halo merger trees finds a roughly constant value 
of (Tl ~ 0.25. In §5, we further discuss the physical significance of cjl and consider whether 
the data can constrain it a priori. Since the Lc is proportional to halo mass, the average 
mass-to-light ratio, T = IO^^'Mq/Liq, is independent of halo mass at a given redshift and 
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uniquely specified by Lio at that redshift. Lio is, therefore, the second free parameter of our 
model. 

For a given halo mass, equations 2.1 and 2.2 give the fraction of halos forming stars 
and the distribution of resulting luminosities, respectively. However, if an unaccounted for 
feedback process further shuts off star formation in some fraction of halos, the resulting 
effective "duty-cycle" will turn our fits of Liq in §4 into underestimates. 

This model gives the average abundance of galaxies in the universe at a given redshift. 
However, there are two geometrical effects that can conspire with structure formation to make 
the observed LF deviate from the average. First, the narrow field-of-view and large halo bias 
at high redshift can result in significant cosmic variance [21-23]. In the z ~ 10 case of Bll, 
the fluctuation in abundance is of order 40%-large but still smaller than other uncertainties in 
the problem. Additionally, the large redshift range probed by high-redshift dropout surveys, 
combined with the rapid evolution of the LF, can have profound implications for the galaxy 
sample, for example, a distribution toward lower redshifts, a higher abundance, and a flatter 
observed LF than would be expected if every galaxy in the sample were at the mean redshift 
[24, see §6] . We show that the contribution of this effect to observations of the z ~ 10 LF 
is much smaller than the uncertainties of the measurements for our model in section §6 and 
ignore it in the rest of this work. 



3 The fitting method 

We fit our model parameters to the data using a minimization method. Focussing on the 
Lio, we marginalize over Msupp, allowing the suppression mass to vary as we search for the 
value of Lio that best fits the data. Thus, our constraints also reflect uncertainties in Msupp. 

Furthermore, we improve upon the MLll constraints by including the upper limits in 
the LFs at z ~ 7 and ~ 8 in our analysis. For each set of model parameters, magnitude 
bins of the LF with measured abundance nobs and error a contribute an amount dx^ = 
[{nmod — nohs) /c]^ to the total value of x^j where rimod is the abundance predicted by the 
model for that set of parameters. Additionally, each magnitude bin in which no galaxies were 
observed and having a la upper limit of n^p contributes dx^ = (?^mod V — 'T-obs / ("-mod V) = 
'T-mod/'T'up) where V is the volume probed by the dropout search in the magnitude bin, and 
we estimate V ^ 1/ n^p . 

At z ~ 10, Bll finds only one galaxy candidate detection (actually an estimate of 
~ 0.8). Therefore, in the associated magnitude bin, a does not include the full Poisson error 
around n-mod for sets of model parameters that produce nmod < n-ohs- Instead, we modify the 
above procedure for this data point so that the error is due to Poisson variance about nmod 
with an additional 50% contribution from cosmic variance. 



4 The evolution of the luminosity function 

In MLll, we presented results for the best-fit values for Msupp at z ~ 6, 7, and 8 and found 
that all were consistent with each other at about 10^'^ Mq. While we were able to constrain 
Msupp strongly from above, a value of 10® Mq could not be ruled out, and the best-fit values 
should be taken as tentative. Here, we provide improved constraints for Lio and extend the 
analysis to z ~ 10. In effect, we extract from the drop in galaxy abundance the contribution 
owing to the rapidly falling HMF and determine whether a complementary change is required 
in Lio, the parameter that describes the average mass-to-light ratio. 
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We consider the observed LF of z ~ 6 i-dropouts [9], z ~ 7 z-dropouts, z ~ 8 Y-dropouts 
[6], and J-dropouts at z ~ 10 [8]. Further improving on the analysis of MLll, we use the 
LFs at z ~ 7 and 8 from the final version of [6], which includes the full two-year HUDF09 
results. We assume a flat prior on Mgupp for values greater than 10^ Mq and fit for Liq 
by minimizing calculated as described in §3. We find values of logio -^^lo = 27.19lQQg^, 
27.36l:[!:J^, 27.69l[]:J^, and 27.45^0^3 at z = 5.9, 6.8, 8.0, and 10.2, respectively. These 
constraints are consistent with the preliminary estimates set in MLll. The rise in Liq from 
z ~ 6-8 represents an average increase in luminosity at fixed halo mass by about 0.24 dex 
per redshift interval. Because of the fitting to the observed LF, the values determined here 
represent observed emission after any extinction by dust. Ref. [9] estimates a ~ 0.18 dex 
dust correction at z ~ 6, implying an intrinsic value of logio -^lo = 27.37 at this redshift. 
On the other hand, the blue spectral slopes of galaxies at z > 7 indicate that they may be 
dust-free [25, but see (author?) 26]. Given the uncertainty, we leave our values uncorrected 
at these higher redshifts. 

The single galaxy candidate observed at z ~ 10 poorly constrains Liq from below. 
However, the lack of detections at brighter magnitudes places interesting limits from above. 
Figure 1 compares the J-dropout LF with the model at z = 10.2 for various combinations 
of Msupp and Liq. The long-dashed curve shows the model LF using the z = 8.0 value of 
Miupp = 10^'^ Mq but extrapolating logio Liq to 28.2 as might be expected if the trend of a 
0.24 dex increase per redshift interval continues to z = 10.2. This model is clearly inconsistent 
with the observations, producing too many galaxies at Muv = —19.6 and perhaps too few 
at Muv = —18.4. Tuning Mgupp cannot reconcile such a large value of Liq with the data. 

On the other hand, the solid and short-dashed lines in Figure 1 assume no evolution in 
the mass-to-light ratio from z = 8.0 and logio Msupp = 9.4 and 8, respectively. While MLll 
showed that the larger value of Mgupp is slightly preferred at z ~ 8, neither is ruled out. Both 
also match the z ~ 10 data well. Given that the model LFs for these values are near the 
two observed upper limits at Muv = —18.9 and -19.6, we conclude that the average galaxy 
brightness for a given halo mass cannot increase as significantly from z ~ 8-10 as between 
z ^ 6-8. 

Figure 2 shows the rise in Liq, corresponding to a drop in the mass-to-light ratio, from 
z ~ 6-10. Because Liq at z ~ 10 is constrained weakly from below, we plot only the upper 
limits. Even correcting for dust at z ~ 6 and ignoring dust at z ~ 8, Liq increases significantly 
between these two redshifts and the respective confidence intervals do not overlap. The 
implied change in galaxy formation over this redshift range contradicts the predictions of 
other studies that the evolution of the LF at z > 6 is entirely due to the behavior of the 
HMF [e.g., 13]. Of course, an individual galaxy does become brighter with decreasing redshift 
[e.g., 15], but its halo mass does not remain fixed. The increase in luminosity with redshift 
at constant halo mass is an important example of how the intrinsic properties of galaxy 
formation are masked by the behavior of the HMF. The population of massive halos shrinks 
faster with increasing redshift than the population of small ones, thus making lower mass 
halos become more typical. This effect "wins" in determining M^y for a Schechter-fitted LF 
even though galaxies at each mass become brighter. The observed decline of M^y is due in 
much larger part to changing halo abundances than to changes in galaxy properties. 

That halos of a fixed mass host brighter galaxies on average at earlier times is expected 
from the increased growth rate of halos and the correspondingly larger galaxy fueling rate 

^ Since the best-fit combination of Lio and Msupp is barely ruled out at the 70% level, the errors here 
represent the 95% confidence interval. 
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Figure 1. The luminosity function of J-dropouts. The point and upper hmits show results from Rcf. 
[8], while the curves show the LF generated with the ML 11 model for different combinations of the free 
parameters Mgupp and Liq. The solid, blue and short-dashed, magenta lines have logio A/gupp = 9.4 
and 8, respectively, with logio ^lo ~ 27.7. The long-dashed, green curve assumes logio (-^supp, -^lo) = 
(9.4, 28.2) with the value of Lio extrapolated to z « 10 based on the 0.24 dex increase per redshift 
interval found from z k, 6-8 (see Fig. 2). The thick solid curve is the same as the thin one except 
that the evolution in the LF along the light-cone (LC) through the Bll redshift selection function 
has been taken into account. 



by cold flows. In §7, we will quantitatively explore the consistency of this prediction. 
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Figure 2. The evolution in the UV luminosity at 10^° from z ^ 6-10. The points show the 
best-fit values at z « 6-8 with the 95% confidence interval (« 2ct) shown at z = 6 and 70% (w Itr) 
shown at z w 7 and 8 (see footnote in main text). The upper magenta point at z « 6 includes a 
dust correction of 0.18 dex. The dashed, green line marks an increase of 0.24 dex per redshift interval 
connecting the points at z w 6 and 8, while the solid, red curves shows the predicted evolution of L\q 
based on the cold-flow model of galaxy fueling and feedback in equation 7.4, assuming a Salpeter IMF 
at solar metallicity with 770 = 1, 4, and 6.3. Blue upper limits denote the 70% and 95% confidence 
constraints placed at z k. 10. 



5 CJl 

The parameter ctl describes the logarithmic spread in the luminosity of galaxies hosted by 
halos of the same mass. In a scenario where the galactic luminosity is due to mergers, ctl 
combines the rate at which galaxies grow fainter by using up gas with the amount of time 
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since the last merger event. For example, if the time since the last merger, tm, in a given 
galaxy is much less than the star formation timescale tgfr; then that galaxy's luminosity will 
be in the bright tail of the distribution in equation 2.2. On the other hand, if tm ^ isfr) the 
galaxy will have used up most of its gas in the starburst and have a luminosity in the faint 
tail of the distribution for its host halo mass. In a cold-flow picture for high-redshift LBGs 
(as further discussed in §7), ctl ~ 0.25 is the result of the factor of a few variation in the 
galactic baryon accretion rate for a give halo mass [11, 27]. 

The value of cl determines how reliably the slope of the LF reproduces that of the HMF. 
For (Tl = with Lc oc Mhaio) every galaxy with the same host halo mass also has the same 
luminosity, so the LF has exactly the same shape as the HMF (at masses sufficiently larger 
than Msupp). Increasing cjl mixes halos of different masses into the same luminosity bin and 
effectively flattens the LF. If we try to fit (Tl a priori from the data (while, for simplicity, 
keeping Mgupp fixed), we find that the constraints from below are poor. This is expected since 
the shape of the LF is not too far from that of the HMF. However, cjl is restricted to be less 
than ~ 0.5 (1), assuming logio-ZVfsupp ~ 9.4 (8). Physically, this implies that the luminosities 
of galaxies hosted in halos of the same mass are all within about an order-of-magnitude of 
each other (or alternatively, within about 2.5 magnitudes of each other. 

6 Observing along the light-cone 

In this section, we discuss the effect of a rapidly evolving LF on a J-dropout survey that 
samples galaxies from a large redshift range. 

The broad photometric bands used in dropout surveys of high-redshift galaxies result 
in significant rerdshift-error. This owes to the uncertain location within the band at which 
a galaxy's spectral break falls. In their Fig. 2, for example, Bll show the redshift selection 
functions of their J-dropout survey spread over a range of almost Az ~ 2. In [24, hereafter 
ML08], we demonstrated that, due to the rapid evolution of the LF, there should be far more 
galaxies in the observed sample at the low-redshift end of the redshift distribution than at the 
high-redshift end simply because the abundance of those galaxies is much higher at a later 
point in the history of the universe. Consequently, the mean redshift of observed galaxies is 
lower than that of the photometric window. For idealized assumptions about the LF and the 
shape of the broad-band filters, we found that the resulting observed abundance of galaxies 
in lO^'^ Mq halos could be about a factor of two higher than the true value expected at the 
mean photometric redshift. 

Here we use the MLll LF model, with logio (Msupp, -Lio) = (9.4,27.7) (see §4), along 
with the exact redshift selection function generated by the Bll J-dropout photometric se- 
lection criteria, to more carefully consider this effect. We calculate the difference between 
the abundance of galaxies at the mean selection redshift of z = 10.2 and that in a survey 
along the line-of-sight light-cone taking into account the evolving LF through the photomet- 
ric band. We find that the density of observed J-dropouts with magnitude Muv = —18.4 is 
about 17% higher than the abundance at z = 10.2 owing to the fact that the mean redshift of 
the observed galaxies is actually z = 9.8 (see Fig. 1). For reference, the star formation rate 
density required to keep the universe ionized is 12% higher at z = 9.8 than at 10.2%. The ef- 
fect increases with increasing brightness causing a slight flatting of the LF (for Muv = —20.4 
the observed abundance is 35% higher than at z = 10.2 and the mean redshift sampled is 
z = 9.7). As discussed in ML08, we ignore the slightly varying differential volume element 
of the survey with redshift as a negligible contribution. 
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While the effect of observing along the light-cone of the Bll survey is not negligible, it 
is less significant than previously anticipated and less than the error already placed on the 
measured abundance. There are three important reasons for this. (1) In ML08, we assumed 
a gaussian distribution for the redshift selection function in ML08 with an exponential tail 
trailing to lower redshifts, while the exact distribution for the Bll J-dropout survey reaches 
down to exactly zero^ at z < 8. (2) The current concordance value of as = 0.82 from the 
Wilkinson Microwave Anisotropy Probe 7-year results [28] is slightly higher than the 3-year 
value of as = 0.776 [29] used in ML08, which results in a HMF that evolves slightly less 
quickly at z ~ 10. (3) The average mass of halos corresponding to a particular brightness 
in the MLll model described in §2 is significantly smaller than in the model of Ref. [12] 
for the LEG mass-to-light ratio that we used in ML08, and so, the abundance of these halos 
varies less rapidly with redshift. Not only is the average brightness of halos at a fixed mass 
larger in the MLll model than in Ref. [12], but the effect becomes even larger when we 
include the decrease in the mass-to-light ratio from z ~ 6, where the model in Ref. [12] was 
calibrated, which was not originally included in ML08. Moreover, the luminosity distribution 
function allows halos of a given mass to dominate the population of galaxies brighter than 
their average; despite the fact that few of these halos shine with luminosities brighter than the 
average, they are much more numerous in total than are more massive halos. These effects 
conspire to slow the evolution of the LF in the range over which the survey is observed. On 
the other hand, if models that predict much higher halo masses for a given luminosity, such 
as the one in Ref. [13], are correct, then the variation in the LF over the redshift range 
of the survey will be substantially larger, and the effects described in ML08 may be quite 
significant. 

7 The evolution of galaxy fueling 

Using a simple but physical model for galaxy formation, we attempt to describe the evolution 
of the average galaxy luminosity for a given halo mass at z > 6. If the star formation time 
scale is much less than the dynamical time on which galaxies grow [30] , then the star formation 
rate is a balance between the rate at which galaxies are fueled by baryons. Mace, and rate 
at which large-scale galactic outflows, ubiquitous at high-redshift [e.g., 31], deplete the gas 
available for star formation. My/. 

MsFR = Mace - Mw. (7.1) 

Here, we have also assumed that a negligible fraction of gas is accreted onto a central black 
hole [32]. 

In the cold flow model, the infall of cold gas onto the galactic disk traces the buildup 
of dark matter in the halo since the gas is never shock-heated to high temperatures [10, 11]. 
Following Ref. [27], the high redshift behavior of this fueling rate is approximated by 

■^The probability that galaxies at 2 < 8 are included in the sample may not be exactly zero but simply 
below the resolution of the Bll simulation depending on the number of mock galaxies used. Since the LF 
increases exponentially with decreasing redshift, even an exponentially vanishing selection probability may 
create a tail of galaxies distributed to much lower redshifts. 
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where /b is the cosmic baryon fraction. Equation 7.2 indicates that gas is accreted more 
quickly at higher redshifts for fixed halo mass. 

According to the momentum-driven wind model of Ref. [33], the wind outflow rate is 
proportional to AfsFR/c, where a is the velocity dispersion. Since a oc M^^^ (1 + z)^/^ [34], 
winds become somewhat less important at higher redshift for fixed halo mass. We thus define 
a mass-loading factor of 

Mw lOOkm/s 
r/w = -r- = m (7.3) 

that describes the relationship between the star formation rate and the wind outflow rate. 
Here, r/o is a free parameter defining the normalization of the wind efficiency with r/o = 1 
in the original derivation from Ref. [33], and ~ 3 a" is the wind velocity at infinity. We 
also note that the recycling of these outflows back onto the galaxy does not affect the baryon 
fueling of equation 7.2 at these early times [35]. 

Assuming a fixed UV luminosity (at 1500 A) to SFR ratio, A = L1500/M* = 8 x 
10^'' erg/s/Hz/(MQ/yr) for a Salpeter IMF and solar metallicity [36, but see also MLll] and 
combining with equations 7.1, 7.2, and 7.3 gives a galaxy UV luminosity as a function of 
mass and redshift for a choice of rjQ. 

T (7 A) 

When ?7o ^ 1 and most of the accreted gas is expelled through winds, L1500 oc M^^^^ (l + z)^, 
while -^1500 Mace when rjQ <^ 1 and most of the accreted gas is turned into stars. In either 
case, the dependence on Mhaio is close to linear, as calculated in the excursion-set analysis 
of MLll. In this model, cjl ~ 0.25 results from the factor of a few spread in the galactic 
accretion rate from galaxy to galaxy around the mean given in equation 7.2 [11, 27]. 

Figure 2 shows the evolution of Liq predicted by the model of galaxy fueling and feed- 
back represented by equation 7.4. Normalizing the model to the value of Liq at z ~ 6 gives 
r/o ~ 4.0 and 6.3 in the dust-corrected and uncorrected cases, respectively. A standard wind 
mass-loading factor with i]q = 1 produces galaxies that are approximately three times too 
bright for their halo mass. This super- wind requirement was previously seen in cosmolog- 
ical numerical simulations [37, 38] which considered r/o = 3 (in our notation). Further, a 
more top-heavy IMF would produce a higher ratio of Lisoo/Mt and require either more dust 
extinction or even stronger winds to reproduce our values of Liq. 

With the model normalized to the value of Lio at z ~ 6, the predicted evolution 
roughly agrees with the behavior we have extracted from the observed LFs. This con- 
sistency is maintained out to z ~ 10, where the a simple extrapolation of the 0.24 dex 
per redshift interval evolution from z ~ 6~8 fails. For rjQ = 4.0 (6.3), the model predicts 
a star formation rate density of 0.03 (0.02) M0/yr/Mpc^ at z = 9.8 (see §6) if M^npp = 
W^Mq and 8(5) x 10""^ Mo/yr/Mpc^ if Msupp = W^-'^ Mq. This compares to a value of 
2.5 X 10^^ C MQ/yr/Mpc^ required to maintain a reionized intergalactic medium [39], 
where /esc is the ionizing escape fraction and C is the clumping factor. 

8 Comparison with recent numerical models 

In this section, we compare our findings with those from several recent numerical and semi- 
analytic models of high-redshift galaxies, all of which produce reasonable agreement with 
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measured LFs. With the flexibihty of our analytic framework, we can cahbrate the physics 
from the observations rather than predict the LF a priori, and our resulting mass-to-light 
ratios consequently "split the difference" between those from other works. 

The numerical simulations of Ref. [15] predict LFs that are systematically brighter 
than observed. While they included some dust extinction calibrated to local measurements, 
excluding dust from their models only exacerbates the discrepancy. Given the blue continuum 
slopes of high-redshift sources [25], assuming additional extinction is probably not justified 
even though it may bring our estimates of the mass-to- light ratio into agreement. The shape 
of their conditional luminosity functions at z = 5.5 are also quite different from ours. While 
our merger-tree analysis predicts a single-peaked distribution (Eq. 2.2), theirs exhibit a 
power-law behavior. At this redshift, our model would indeed anticipate a rising abundance of 
fainter galaxies hosted by 1.7 x 10-*^° Mq and 3.4 x W^^ Mq halos at least down to Mjjy = -18 
(i.e., the faint limit of their Fig. 14). However, we would have expected the distribution for 
1.6 X 10^^ Mq halos to have turned over by around Muv = —19.5. The contribution from 
faint satellites, neglected in our treatment, may be responsible. We note that their simulation 
includes momentum-driven winds with mass-loading factors similar to what we found in §7. 

On the other hand, the simulations by Ref. [14] generate galaxies at z = 7 that are 
about a factor of two (i.e., about 0.7 magnitudes) more luminous than ours for 10^^ 
halos. Yet, curiously, their resulting LF at this redshift is shifted somewhat from the data 
toward fainter magnitudes. Some of this discrepancy may result from their supernova wind 
prescription, which assumes an initial wind velocity of 500km/s for all halos even though 
supernovae may not expand very far into the dense interstellar media of high-redshift galaxies 
[40]. It is difficult, however, to compare their results either to ours or to observations given 
their relatively small box size (only 10/i^^Mpc), which contains a very limited sample of 
simulated bright galaxies. Still, their detailed treatment of metal-enrichment may make 
their conclusion of a negligible dust contribution to the LF fairly reliable. 

Finally, Ref. [41] use a semi-analytic model of galaxy formation, combining N-body 
simulations of dark matter with analytic prescriptions for star formation, to model the high- 
redshift LF. Their results are significantly different from both Ref. [15] and Ref. [14] in 
that both mergers and a topheavy IMF are required to reproduce observations. However, it 
is possible that their need for a topheavy IMF is a result of over-quenching by their strong 
feedback prescription from supernovae-driven winds. At the same time, their requirement 
for mergers may result from the long timescale for quiescent star formation assumed from 
low-redshift observations; if high redshift galaxies have higher molecular fractions for the 
same total gas mass [e.g., 42], star formation will proceed more efficiently that this estimate 
[e.g., 43]. 

9 Discussion and conclusions 

We have investigated the observed evolution in the LF of LBGs from z ~ 6~10 in the 
context of a rapidly changing halo mass function to disentangle the corresponding variation 
in galaxy formation. Our analysis improved on our previous work in MLll by considering 
LFs based on the complete two-year HUDF09 data as well as incorporating information from 
null-detections at bright magnitudes. 

From z ~ 6-8, we find a significant rise in LiQ-the luminosity at 1500 Aof a galaxy 
hosted in a halo of lO^'^ Mq and the parameter that describes the average galaxy mass- 
to-light ratio-which is inconsistent with the predictions of other analytic models [e.g., 13]. 
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Despite the detection of only a single J-dropout at 2; w 10, we were also able to place 
interesting constraints on galaxy formation at this high redshift. The upper limits on the 
number of galaxies with Muv ^ ~19 were particularly useful in this regard. We were able to 
rule out a decrease in the average mass-to-light ratio of galaxies from z ~ 8-10 comparable 
to the evolution from z ~ 6-8. This result would remain unchanged if the single Bll galaxy 
candidate were actually a spurious detection or a low-redshift interloper (Bll estimate a 20% 
total probability for such an event). However, there are two observational effects that could 
allow for a higher halo luminosity (i.e. lower mass-to-light ratio) at z ~ 10. First, the strict 
set of selection criteria designed to minimize low-redshift interlopers could have resulted in 
undetected sources that should be considered in the LF [see, e.g., 44]. Second, light outside 
the selection aperture could have been missed due to surface brightness dimming; the Bll 
object and upper limits were not corrected for this potential missing signal. The correction 
can increase the brightness by ~ 0.5 magnitudes at z ~ 8 [6]. Since the surface brightness 
of a galaxy with fixed size and spectral luminosity decreases"^ as (1 + z)~^ and assuming 
galaxy size scales as (1 -|- 2;)~o.875 j^-j^j^ estimate that the correction at z ~ 10 could be 
as large as a full magnitude making the data not inconsistent with the 0.24 dex per redshift 
trend shown in Figure 2. Additionally, if we have failed to account for a feedback process 
that shuts off star formation and creates an "effective duty-cycle," reducing the abundance of 
halos seen as galaxies, then the true values of Lio will be higher than we have calculated with 
the degree of underestimation depending on the mass, luminosity, and redshift dependencies 
of the feedback mechanism. 

While we calculated the above results with ctl = 0.25 based on the MLll merger tree 
calculation, we also examined to what extent this value is independently constrained by the 
data if left to vary. An increase in ctl flattens the slope of the LF by mixing halos from a 
larger mass range into the same luminosity bin. We find that the data limit the logarithmic 
spread in galactic luminosity at fixed halo mass to be less than about an order-of-magnitude. 

We then presented a simple but physical model for galaxy fueling from cold-flows and 
feedback from momentum-driven winds to explain the evolution we found in L\q. Normalizing 
to Lio at z ~ 6, this model does an excellent job of predicting the behavior of the average 
galaxy luminosity out to 2; ~ 10. Thus, it appears that, on average, galaxy formation at these 
redshifts can be described as a steady-state balance among gas accretion, star formation, 
and wind ejection, as discussed in [30], while ignoring the growth of the gas reservoir in the 
disk. Moreover, if galaxy formation proceeds in halos as small as 10* Mq [34], then the star 
formation rate density at z ~ 10 should be high enough to keep the universe ionized as long 
as/e^c'C< 10. 

However, future data from WFC3, both deeper and in multiple fields (to reduce cosmic 
and Poisson variance), will improve constraints on the observed LF and further test these 
predictions even before the James Webb Space Telescope comes online. 
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